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ABSTRACT 

We report on IRAC-4.5 /Ltm, IRAC-8.0 /im and MIPS-24/im deep observations of 16 Gamma-Ray 
Burst (GRBs) host galaxies performed with the Spitzer Space Telescope, and we investigate in the 
thermal infrared the presence of evolved stellar populations and dust-enshrouded star-forming activity 
associated with these objects. Our sample is derived from GRBs that were identified with sub- 
arcsec localization between 1997 and 2001, and only a very small fraction (~20%) of the targeted 
sources is detected down to /4.5^ m ~ 3.5 /xJy and f24p,m ~ 85 /iJy (3cr). This likely argues against a 
population dominated by massive and strongly-starbursting (i.e., SFR > 100 M Q yr -1 ) galaxies as it 
has been recently suggested from submillimctcr/radio and optical studies of similarly-selected GRB 
hosts. Furthermore we find evidence that some GRBs do not occur in the most infrared-luminous 
regions - hence the most actively star- forming environments - of their host galaxies. Should the 
GRB hosts be representative of all star-forming galaxies at high redshift, models of infrared galaxy 
evolution indicate that > 50% of GRB hosts should have /24/jm <; 100 /iJy. Unless the identification of 
GRBs prior to 2001 was prone to strong selection effects biasing our sample against dusty galaxies, 
we infer in this context that the GRBs identified with the current techniques can not be directly used 
as unbiased probes of the global and integrated star formation history of the Universe. 

Subject headings: galaxies: high-redshift — infrared: galaxies — cosmology: observations — galax- 
ies: individual (GRB970828, GRB980425, GRB980613, GRB980703, GRB981226, 
GRB990705) 



1. INTRODUCTION 

It is now widely believed that the so-called "long" 
Gamma-Ray Bursts (i.e., GRBs with duration >2s and 
soft spectra, as opposed to the short and hard bursts, 
e.g., Kouveliotou et al. 1993) are intimately connected 
to the collapse and the cataclismic destruction of some 
short-lived and very massive stars. Evidence in this re- 
gard include (i) the signature of Type Ic supernova and 
hypernova in the optical transient emission of GRB coun- 
terparts (e.g., Galama et al. 1998; Stanek et al. 2003; 
Hjorth et al. 2003; Malesani et al. 2004); (ii) heavy el- 
ements from metal-enriched media typical of supernova 
remnants observed in the spectrum of several X-ray af- 
terglows (Piro et al. 2000; Reeves et al. 2002); (iii) the 
unquestionable starbursting nature of their host galax- 
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ies (e.g., Bloom et al. 1998; Fruchter et al. 1999; Sokolov 
et al. 2001; Chary et al. 2002; Le Floc'h et al. 2003; Chris- 
tensen et al. 2004); and (iv) the location of GRBs relative 
to the center of their hosts, which appears to be consis- 
tent with a population of progenitors residing in galaxy 
disks (Bloom et al. 2002). Hence it has often been pro- 
posed that long GRBs could be used as powerful tracers 
of the global star-forming activity in the early Universe 
(e.g., Wijers et al. 1998; Mirabel et al. 2000; Blain & 
Natarajan 2000). As illustrated by the recent identifica- 
tion of a burst at z = 6.29 (e.g., Kawai et al. 2005), GRBs 
are indeed detectable up to very high redshifts (Lamb 
& Reichart 2000). They are also very little affected by 
dust extinction, which is known to be particularly signif- 
icant in distant starburst galaxies (e.g., Blain et al. 1999; 
Franceschini et al. 2001; Chary & Elbaz 2001). 

However, this picture relies on the strong assumption 
that the production of GRBs in a given starburst region 
scales only with the rate of stars formed in this environ- 
ment, with no dependency on other physical parameters 
that may vary from one galaxy to another or that may 
evolve throughout the lifetime of the Universe. This as- 
sumption is rather major, as the occurrence of a GRB 
could strongly depend on the properties of the gas from 
which its progenitor originates (e.g., metallicity: Mac- 
Fadyen & Woosley 1999; Ramirez-Ruiz et al. 2002; Heger 
et al. 2003; Hirschi et al. 2005). It could also vary accord- 
ing to the fraction of those progenitors involved in binary 
systems (Izzard et al. 2004; Podsiadlowski et al. 2004; 
Mirabel 2004b), and it finally relies on a non-evolution 
• r 'of the Initial Mass Function with redshift. Testing this 
"one-to-one" connection between GRBs and star forma- 
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tion is therefore particularly crucial to guarantee an ac- 
curate understanding of the use of GRBs as quantitative 
tracers of galaxy evolution. 

One possible approach to investigate this relation is 
to compare the properties of the hosts of Gamma-Ray 
Bursts with respect to the galaxies responsible for the 
bulk of the star-forming activity in the Universe as a 
function of redshift. It is now well established that a sig- 
nificant fraction of the present-day stellar mass budget 
was formed during brief and violent infrared-luminous 
(Lir = L 8 _ioooMm 10 11 L Q ) episodes of star formation 
within massive (M >5x 10 10 M ) galaxies at 0.5 < z < 3 
(e.g., Flores et al. 1999; Blain et al. 2002; Elbaz et al. 
2002; Dickinson et al. 2003; Lagache et al. 2003; Frances- 
chini et al. 2003; Le Floc'h et al. 2005; Hammer et al. 
2005). Investigating the extinction-corrected star forma- 
tion rate and the mass of GRB hosts should therefore 
provide tight constraints on the relevance of GRBs for 
probing the star formation history of the Universe. 

Previous observational studies on GRB hosts though 
have led to conflicting views about their nature. Based 
on their properties at optical and near-infrared wave- 
lengths it has been argued that GRB hosts are mostly 
blue, sub-luminous and low-mass galaxies with young 
stellar populations, characterized by a modest activity 
of star formation and potential selection effects due to 
low metallicity (e.g., Sokolov et al. 2001; Le Floc'h et al. 
2003; Fynbo et al. 2003; Courty et al. 2004; Prochaska 
et al. 2004; Christensen et al. 2004). On the other hand 
it has also been claimed that the morphology of these 
objects and their average radio/submillimcter properties 
rather indicate massive and actively star-forming galax- 
ies (Conselice et al. 2005; Berger et al. 2003). In order to 
bring tighter constraints on the nature of those sources, 
we undertook a survey of GRB host galaxies with the 
Spitzer Space Telescope. In this paper we present mid- 
infrared (mid-IR) images of 16 objects at 4.5 /zm, 8 /xm 
and 24 ^m, which allows us to constrain the presence 
of evolved stellar populations and dust-enshrouded star- 
forming activity in these sources. In Sect. 2 we describe 
the Spitzer data used in this study. Sect. 3 outlines 
general results derived from these mid-IR observations, 
while the GRB host properties are detailed on a galaxy- 
by-galaxy basis in Sect. 4. Our findings are discussed 
in Sect. 5, and they are finally summarized in Sect. 6. 
Throughout this work, we assume a ACDM cosmology 
with H = 70 km s^Mpc" 1 , ft m = 0.3 and Q x = 0.7 
(Spergel et al. 2003). 

2. OBSERVATIONS 

2.1. The data 

With the exception of the GRB 970228 which was de- 
tected in a sub-luminous blue dwarf galaxy at z = 0.695 
(Bloom et al. 2001a), our sample of GRB host galax- 
ies was built by considering every GRB identified before 
July 10 th , 1999 and localized on the sky with a sub- 
arscsec accuracy thanks to the detection of an afterglow 
at longer wavelengths. This selection is therefore inde- 
pendent of any a priori information on the properties 
of the hosts (e.g., redshifts, star- formation rates, lumi- 
nosities, ...). Furthermore it does not take into account 
at which wavelength the position of the GRB was deter- 
mined, i.e., whether the afterglow was optically bright or 



whether it was only seen in the radio and/or the X-rays 
(i.e., "dark" burst). This led to a sample of 15 objects, to 
which we added the host of the burst detected on Febru- 
ary 22 nd , 2001 (i.e., GRB 010222). This host galaxy has 
been claimed to be associated with a SCUBA/MAMBO 
ultra- luminous infrared object (Frail et al. 2002), making 
it an obvious and quite interesting target for IR observa- 
tions. 

We note however that our sample is obviously subject 
to the various observational cuts that may have biased 
the identification of GRBs prior to 2001. For instance, 
most of the bursts considered in this work were de- 
tected either with the WCF camera on-board the Bcppo- 
SAX satellite or with the BATSE instrument on-board 
the Comptom Gamma-Ray Observatory. As revealed 
by more recent high energy missions like HETE-2 and 
Swift , this may have imposed a first sub-selection over 
the whole population of GRBs depending on their intrin- 
sic luminosity and the hardness of their spectrum. Fur- 
thermore, the subsequent ground-based follow-ups that 
led to the sub-arcsec localization of these GRBs as well 
as the determination of their redshift from optical spec- 
troscopy were likely slower than those currently oper- 
ated by the new networks of more dedicated telescopes. 
This could in principle favor the detection of X-ray and 
optical transients relatively brighter than the afterglows 
accompagnying the typical bursts now accessible with 
Swift (Berger et al. 2005), thus biasing the selection to- 
ward host galaxies with low extinction. In Sect. 5.4 we 
will discuss in more detail the implication of our results 
taking into account the above mentioned issues. 

The observations were performed with the Spitzer 
Space Telescope (Werner et al. 2004) as part of the 
IRS Guaranteed Time Observing program (Houck et al. 
2004). The J2000 coordinates and the redshifts of the 
targeted galaxies are reported in Table 1 along with the 
name of the GRBs that led to their selection. We ini- 
tiated our Spitzer program in 2004 while these gamma- 
ray bursts all occured before February 2001. The con- 
tribution of the emission from their fading afterglow was 
therefore negligible at the time of our observations. 

Each object was imaged with the "InfraRcd Array 
Camera" (IRAC, Fazio et al. 2004) at 4.5 and 8.0 /j,m 
as well as with the "Multi-band Imager and Photometer 
for Spitzer" (MIPS, Rieke et al. 2004) at 24 ^m. The 
IRAC detectors are characterized by 256x256 squared 
pixel arrays with a pixel size of 1.22" leading to a total 
field of view of 5.2x5.2 arcmin 2 . The Full- Width at Half 
Maximum (FWHM) of the Point Spread Function (PSF) 
varies between 1.5" and 2" accross the different IRAC 
channels. The MIPS detector at 24fim uses a 2.45" pixel 
size array of 128 x 128 elements also resulting in a field of 
view of 5.2x5.2arcmin 2 . The image at this wavelength 
is characterized by a PSF with a FWHM of ~6". 

Each observation was performed using a sequence of 
several frames slightly dithered with respect to the po- 
sition of the GRB host. Three frames of 100 s were ob- 
tained at 4.5 /im and 8 /im giving a total integration time 
of 300 s per source and per band. At 24 /im, 14 frames of 
30 s led to a total exposure time of 420 s per source. The 
data were reduced with standard procedures (i.e., dark- 
current subtraction, cosmic ray removal, non-linearity 
correction, flat-fielding and mosaiking) using the pipeline 
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of Spitzer Science Center 8 . 

The absolute pointing accuracy of the Spitzer satellite 
is better than ~ 1". The la relative astrometric uncer- 
tainty is less than - 0.3" in the IRAC and MIPS data. 

2.2. IRAC and MIPS Photometry 

With the exception of the GRB 980425 host galaxy 
lying at z = 0.0085 (Galama et al. 1998; Tinney et al. 
1998), the other GRB hosts of our sample are all lo- 
cated at cosmological distances (z>0.84, see Table 1). 
Therefore they are not spatially resolved in the Spitzer 
images and their fluxes can be estimated using small cir- 
cular aperture photometry. In the IRAC images, counts 
were measured over a circled area with a radius of 3 pixels 
(e.g., 3.6") centered at the position of each target. These 
counts were translated into flux densities using the con- 
version factor prescribed in the Spitzer Observing Man- 
ual 9 and a slight correction was finally applied to account 
for the extended size of the PSF. Sensitivity limits of 
~3.5 /iJy and ~ 20 /iJy (3a) were estimated at 4.5 /im and 
8.0 fim respectively, based on the dispersion of the flux 
measurements obtained over blank sky regions within the 
same aperture as the one used for the photometry of the 
objects. 

Evidence for non-negligible blending was found at 
24 ^m, which is mostly due to the larger size of the 
PSF and the higher level of extragalactic confusion in 
the MIPS images than in the IRAC data. To ensure ac- 
curate results, source extraction and photometry were 
therefore performed using the PSF fitting technique of 
the DAOPHOT software (Stetson 1987). We constructed 
an empirical point spread function from the brightest 
point sources found in our 24 /im data. This PSF was 
accordingly scaled to provide the best match to the ob- 
ject detected at the position of the GRB host, thus lead- 
ing to a direct estimate of its total flux at 24 ^im. As 
a sanity check, the same procedure was also performed 
using a 24^m theoretical PSF simulated and provided 
by the Spitzer Science Center. Within the uncertainties, 
the photometric measurements that we obtained in this 
case are consistent with those derived using the empirical 
point-spread function. A 3a sensitivity limit of ~ 85 fiJy 
was estimated using aperture photometry 

Regarding the nearby GRB 980425 host galaxy, a dif- 
fuse extended infrared emission was detected up to a dis- 
tance of 20" to 30" from the center of the object. The 
sky background level was therefore estimated within an 
annulus defined by an inner radius of 35" and a width 
of 10". At 4.5 ^m the total flux of the galaxy was deter- 
mined using a 23"-radius aperture. This allowed us to 
recover most of the extended emission of the object while 
avoiding the contamination from other field sources lo- 
cated close to the host. At 8 /tm and 24 ^m, a 35"-radius 
aperture was found to provide a very good estimate of 
the total flux of the galaxy. At these wavelengths the 
source density in the field is smaller than observed at 
4.5 ^m, and no other contaminant object was detected 
within this large aperture. 

Our flux measurements and upper limits are given in 
Table 1. The absolute photometric uncertainties in the 

8 see http://ssc.spitzer.caltcch.edu/postbcd/ 

9 An electronic version of the Spitzer Observing Manual is avail- 
able at http://ssc.spitzer.caltech.edu/documcnts/SOM/ 



IRAC and the MIPS data are respectively less than 5% 
and 10%. 

3. RESULTS 

3.1. Detection rate of the GRB host galaxies in the 
IRAC and MIPS images 

The 4.5/8.0 /xm IRAC and 24 ^m MIPS images of the 
GRB host galaxies for which a detection in at least one 
mid-IR band was obtained (see also Table 1) are pre- 
sented as postage stamps in Figures 1 to 6. To facilitate 
the identification of the fields of view, they are displayed 
along with optical images publicly available in the liter- 
ature. These optical data were obtained with the STIS 
or the WFPC cameras on-board the Hubble Space Tele- 
scope (HST) as part of various observing programs led 
by Fruchter et al. (HST Proposals 7966/8189), Holland 
et al. (HST Prop.: 8640) and Kulkarni et al. (HST 
Prop.: 8867). Most of the reduced images are taken 
from the "Survey of the Host Galaxies of Gamma-Ray 
Bursts" 10 (Holland et al. 2000a), with the exception of 
the data for the hosts of GRB 970828 and GRB 010222. 
For the latter, reduced products were provided by the 
"Multimission Archive at the Space Telescope Science 
Institute" (MAST) 11 . Most of these HST images are dis- 
played with higher spatial resolution in the Figure 2 of 
Bloom et al. (2002). 

The nearby host of GRB 980425 is clearly detected 
in the three bands that we covered with IRAC and 
MIPS. However, we note that most of the cosmologi- 
cal GRB hosts from our sample are not detected with 
Spitzer. Among the 15 high-rcdshift targeted sources, 5 
are brighter than 3a at 4.5 /im, while only 1 and 3 ob- 
jects are detected at 8/Ltm and 24 /im respectively. At 
z ~ 1 a typical L* galaxy is most often an intermediate- 
mass IR- luminous spiral (e.g., Zheng et al. 2004; Ham- 
mer et al. 2005; Le Floc'h et al. 2005; Melbourne et al. 
2005), with an SED leading to observed flux densities of 
f 4. 5 am ~ 10 [i Jy, /spm ~ 8 /«Jy and f 2 A^ m ~ 90 ^tJy. Given 
the 3a sensitivity limits of our IRAC and MIPS data (see 
Sect. 2.2), these non-detections consequently have strong 
implications on the nature of the GRB hosts relative to 
typical field star-forming galaxies. In Sect. 4 we will de- 
scribe on a case-by-case basis the constraints that can be 
derived from these results regarding the spectral energy 
distributions (SEDs) of each galaxy. We will discuss their 
global implications for our general understanding of the 
GRB hosts in Sect. 5. 

For each of the three Spitzer bands, the non-detected 
sources were also stacked together in the attempt of in- 
fcring deeper constraints on the average mid-IR fluxes of 
these objects considered as a whole population. However 
the small number of stacked images improved the original 
depth of our data by only a factor of ~ 1.5-3 depending 
on the wavelength, and no signal was detected above the 
resulting 3a levels. We note that even though we per- 
formed the stacking test, this approach should ideally be 
applied only for sources located in a thin redshift slice, so 
that the constraint on the flux measured in the stacking 
can be converted into a more physical quantity associated 
with these objects. As a result, the broad range of high 
rcdshifts covered by our targets (see Table 1) prevents a 

10 http:/ /www. ifa.au.dk/~hst/grb_hosts/intro. html 

11 http://archive.stsci.edu 
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robust interpretation of our lack of detection in the fi- 
nal stacked GRB host images. It does suggest though 
that considerably deeper observations will be required in 
order to increase the rate of GRB host detections with 
Spitzer. 

3.2. Infrared luminosities 

In the local Universe, correlations between the mid- 
IR luminosity and the 8-1000 fim integrated emission of 
galaxies have been observed not only for normal and qui- 
escently star-forming objects (Dale et al. 2001; Rous- 
scl et al. 2001) but also for more actively starbursting 
sources (Chary & Elbaz 2001). Such correlations likely 
hold also in the distant Universe, as the mid-IR/radio 
relationship resulting from the mid-IR/far-IR and far- 
IR/radio local correlations is still observed at high red- 
shifts (Gruppioni et al. 2003; Appleton et al. 2004). 
Therefore the mid-IR Spitzer data can potentially be 
used to infer some constraints on the total infrared lumi- 
nosity of distant galaxies. 

In this extrapolation, the errors are largely dominated 
by the uncertainties on the shape of the underlying SED 
from 8 ^(m to 1000 ^m. For example, the predictions from 
the libraries of starburst-dominated SEDs proposed by 
Chary & Elbaz (2001), Dale et al. (2001) and Lagache 
et al. (2003) appear to be consistent within only 0.3 dex 
up to z ~ 2, but this dispersion can be significantly larger 
assuming other SEDs (e.g., Dale et al. 2005). At z> 1.5 
in fact, our 24 /tm data are only sensitive to the brightest 
galaxies (i.e., ULIRGs), which harbor a large diversity of 
mid-IR properties (e.g., Armus et al. 2004, Armus et al. 
2006 in prep.). The prediction of their total IR luminos- 
ity based on a single mid-infrared flux measurement can 
be uncertain by a factor of > 5 for a given object. 

To constrain the bolometric luminosity of the GRB 
host galaxies, we first converted the 24 /<m flux density 
(or upper limit) measured for each object with a known 
spectroscopic redshift to a monochromatic luminosity at 
the corresponding rest-frame wavelength 24/xm/(l+z). 
Following the approach presented by Le Floc'h et al. 
(2005), this estimate was then translated into a total 
IR luminosity using the three libraries mentioned in the 
previous paragraph. We should emphasize that in these 
collections of IR spectra, a given monochromatic lumi- 
nosity at a given wavelength corresponds to a single "to- 
tal IR" luminosity, thus leading to a unique determina- 
tion of this quantity for each 24 /jm flux measured in our 
sample. The estimates derived from the three libraries 
were therefore averaged for each object and their dis- 
persion was used to quantify the associated uncertainty. 
We consider our results to be accurate within a factor of 
~ 2-3 up to z ~ 1.5 and within a factor of ~ 3-5 at higher 
rcdshifts. We also keep in mind that such uncertainties 
apply on a case-by-case basis, and they are obviously 
smaller when addressing the average IR luminosity of a 
sample of galaxies. 

These estimates are reported in Table 1 along with the 
equivalent star formation rates (SFR) derived from these 
infrared luminosities using the calibration proposed by 
Kennicutt (1998). These conversions assume that the 
whole IR emission detected at 24 /im originates from star- 
forming activity. While the host of GRB 980425 is only 
a modest infrared emitter (Lir = 2 x 10 9 L Q ), the three 
24 /mi-detected high- redshift sources (i.e., the hosts of 



GRB 970828, GRB 980613 and GRB 990705) are charac- 
terized by a total infrared luminosity Lir in the range 
of 10 11 L Q < Lir < 10 12 L , bringing them to the class 
of the so-called Luminous InfraRed Galaxies (LIRGs, 
Sanders & Mirabel 1996). However, there is no detec- 
tion of a host with luminosity larger than 10 12 L© in the 
range of the ULIRGs (Ultra-Luminous Infrared Galax- 
ies: 10 12 L < Lir < 10 13 L ) and the HyLIRGs (Hyper- 
Luminous Infrared Galaxies: Lir > 10 13 L©). In the case 
of the non-detections, the measured 3a sensitivities were 
used to derive an upper limit on the total IR luminosity 
when a confirmed redshift was available. Most of these 
constraints are also consistent with infrared luminosities 
lower than 10 12 L©, even though we can not definitely 
rule out having a few GRB hosts with ULIRG-typc IR lu- 
minosities given the uncertainties affecting our estimates. 

One might note a few potential caveats that could af- 
fect these total IR luminosity estimates. First, we can 
not completely exclude the presence of a dust-embedded 
active galactic nucleus (AGN) lurking in these GRB host 
galaxies and dominating their mid-IR emission. Because 
the IR SED of AGNs is generally much flatter than the 
starburst-dominated SEDs that we assumed in the con- 
version of the 24 fim flux density (see Weedman et al. 
2005), our determination of the total IR luminosity could 
be over-estimated by a factor of 5 to 10 in these cases. 
However, we consider this possibility unlikely as no typ- 
ical AGN signature has ever been reported from the op- 
tical, X-ray and radio properties of these objects. 

Furthermore, we have implicitly assumed that the 
MIPS 24 ^im detections purely originate from star- 
forming activity in the hosts and we have neglected a 
possible contribution from a transient emission due to 
the effect of the GRB on its close environment. Based 
on a detailed modeling of the heating effect of GRBs 
and their afterflows on their surrounding region, Vene- 
mans & Blain (2001) have argued that a reprocessed dust 
mid-IR emission could in principle be easily detected by 
Spitzer several years after a burst occuring in a dusty 
star-forming galaxy. In our analysis of the IR spectral 
energy distributions of the GRB hosts, presented in Sec- 
tions 4 & 5, we do not consider this likelihood. We will 
argue though that the contribution of such a GRB dust 
emission is unlikely to be significant in our sample. 

3.3. Spectral energy distributions 

The GRB host galaxies of our sample have already 
been extensively observed at optical and near-infrared 
wavelengths, and some of them have also been targeted 
by submillimeter and radio observations. Consequently, 
our Spitzer data can be used in a multi-wavelength con- 
text to derive some constraints on the global spectral 
energy distribution and the nature of these sources. 

We combined our mid-infrared photometry with other 
broad-band imaging data retrieved from the literature 
(see the caption of Fig. 7 for references). Optical and 
near-infrared magnitudes were first corrected from the 
foreground Galactic extinction using the DIRBE/IRAS 
dust maps of Schlegel et al. (1998) and assuming the 
iiy = 3.1 extinction curve of Cardelli et al. (1989). De- 
pending on the quoted magnitude reference, they were 
then converted into fluxes using the zero points from the 
Vega or the AB systems. These fluxes as well as the 
fluxes or the upper limits gathered at the other wave- 
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lengths were finally converted into rest-frame monochro- 
matic luminosities to provide constraints on the GRB 
host SEDs from the optical up to the mid-IR or the sub- 
millimeter/radio wavelength range. 

Our results are illustrated in Figure 7. With the ex- 
ception of the GRB 990510 host galaxy that we did not 
detect with Spitzer and which is also optically very 
faint (V~ 28 mag, Fruchtcr ct al. 2000a), all GRB hosts 
with a confirmed spectroscopic redshift in our sample 
are shown. Given the very low rate of detections at 
long wavelengths (i.e., A>3^m), we do not perform 
a fit to the current data and we refer the reader to 
e.g., Sokolov et al. (2001), Christensen et al. (2004) and 
Chary et al. (2002) for a statistical stellar population 
synthesis of the GRB host optical and near-IR proper- 
ties. Rather, we consider as a comparison the spectral 
energy distribution of the prototypical objects NGC 253, 
M82 and Mrk231. NGC 253 and M82 are moderately- 
active star-forming spirals characterized by IR luminosi- 
ties Lir ~ 1.4x 10 10 Lq and L m ~ 3.7x 10 10 L Q respec- 
tively (Forster Schreiber et al. 2003). Their mid-IR SEDs 
were derived from the ISOCAM-CVF spectra of Forster 
Schreiber et al. (2003), which were further extrapolated 
to the UV/optical and the far-IR/radio domains using 
the modelling provided by Silva et al. (1998). Mrk231, 
on the other hand, is a warm ULIRG with an IR lu- 
minosity Lif{/^ 4x 10 12 Lq powered both by a luminous 
AGN and a violent starburst (Farrah et al. 2003). Its 
SED was derived by combining the 5-35 /zm spectrum 
recently obtained by the IRS on-board Spitzer (Weed- 
man et al. 2005; Armus et al., in prep.) with far-IR 
and radio data published in the literature (e.g., Ivison 
et al. 2004). In addition to these three sources, we also 
displayed the template of a starburst-dominated ULIRG 
with Lir = 4 x 10 12 Lq as well as the SED of a cold LIRG 
with Lir, = 10 11 L . These two spectral energy distribu- 
tions were taken from the IR galaxy libraries derived by 
Chary & Elbaz (2001) and Lagache et al. (2003) respec- 
tively. 

All these SEDs were chosen to be globally representa- 
tive of the expected emission from starburst galaxies and 
IR-luminous sources of the local Universe. They are dis- 
played in Figure 7 with the y-axis in units of WHz -1 in 
order to provide a simple and direct qualitative compar- 
ison with the observed properties of the GRB hosts as a 
function of wavelength. As a result they should only be 
viewed in this context. Since no fitting to the detections 
and upper limits was attempted, some of these SEDs 
were actually re-scaled to match the luminosity of the 
GRB hosts at certain wavelengths. This normalization 
was done in an ad-hoc manner, either at the shortest 
detected wavelengths (e.g., B-band), in the rest-frame 
near-IR (e.g., if-band or IRAC 4.5 ^.m channel), at the 
MIPS 24 /im observed wavelength, or even in the radio 
as in the case of the GRB 980703 host. It better re- 
veals how the global SEDs of the GRB-selected galaxies 
deviate from the other templates considered in the fig- 
ure, and it also shows the large diversity of properties 
characterizing the GRB host population. This will be 
more thoroughly discussed on a galaxy-by-galaxy basis 
in Sect. 4. 

4. THE GRB HOSTS IN A MULTIWAVELENGTH CONTEXT 



Following the results derived in Sect. 3.3 we analyze 
hereafter the multi-wavelength properties of the GRB 
host galaxies using the spectral energy distributions dis- 
played in Fig. 7. We first discuss on a source-by-source 
basis the objects with a known spectroscopic redshift and 
ordered by increasing distance from Earth, and we briefly 
mention the properties of the remaining sources at the 
end of the section. 

4.1. The host galaxy of the GRB 980425 

The association of the GRB 980425 with the 
SN1998bw initially proposed by Galama et al. (1998) 
led to the identification of a nearby sub-luminous blue 
galaxy located at z=0.0085. Classified as an Sbc type in 
the Hubble sequence (Fynbo et al. 2000), it is currently 
the only GRB-selected object known in the local Uni- 
verse. The supernova SN 1998bw was observed in one of 
its spiral arms at a distance of ~ 900 pc from a rather 
bright HII region. 

As shown in Fig. 2 this galaxy is detected with a high 
signal to noise in the three Spitzer bands. Because of 
its proximity it is also well resolved and a diffuse emis- 
sion extending along the major axis is clearly observed. 
Nonetheless the most striking result is the detection of 
a bright mid-IR point source located very close to the 
region where the GRB occured. It is the brightest point 
source detected at 4.5 /im and its contribution to the to- 
tal emission of the galaxy is steeply rising with wave- 
length. With a flux density of 220 /iJy and 1815 /xJy at 
4.5 /xm and 8/zm, it represents respectively 7% and 15% 
of the monochromatic luminosity of the whole galaxy 
measured in the IRAC channels. At 24 fim its flux 
reaches ~21mJy, implying that more than ^75% of the 
energy radiated by the galaxy at the MIPS wavelength 
arises from this region. 

This luminous point source represents one the reddest 
objects so far identified in the nearby Universe. With a 
24 ^m monochromatic luminosity ^24^111^1-1 x 10 8 Lq, 
it is also much more luminous than W49 which is the 
brightest HII region within the Milky Way (Harper & 
Low 1971). As a result, its discovery in a GRB-selected 
galaxy raises obviously the question of a physical link 
with the hypernova. In spite of the large PSF character- 
izing the Spitzer data, the brightness of this object al- 
lowed us to estimate its position along the spiral arm of 
the galaxy with an accuracy better than 0.4". We believe 
that it is unlikely to be related to the close environment 
of the GRB but it rather coincides with the bright HII 
region located ~5" (i.e., ~900pc) in the North- West di- 
rection. It could be due, for instance, to a dense super 
star cluster deeply embedded in dust. This source is to 
be observed very soon with the Infrared Spectrograph 
of Spitzer to explore in more detail its nature and its 
mid-IR properties, and the results will be reported in a 
forthcoming paper. 

Should this object be indeed related to intense dusty 
star formation, our mid-IR observations as well as the 
optical view of the galaxy reveal that the GRB did not 
occur in the most active site of star formation within its 
host. 

4.2. The host galaxy of GRB 970508 

This object is a blue compact dwarf galaxy (Fruchtcr 
et al. 2000b) located at z = 0.835 and showing at op- 
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tical wavelengths a quiescent activity of star formation 
(Bloom et al. 1998). Its 4000-8000A spectrum reveals a 
blue continuum with an [Oil] emission line corresponding 
to a star formation rate SFR [ OII ] ~ 1M Q yr _1 (not cor- 
rected for dust extinction, Bloom et al. 1998). It is not 
detected with IRAC nor with MIPS. Given its moderate 
rcdshift and the depth of our data, this non-detection 
with Spitzer and the flux density measured in the K- 
band exclude the presence of a massive underlying stel- 
lar population in this object. It also argues against a 
significant contribution from dust-obscured star-forming 
activity. As a result the total star formation rate is un- 
likely to be much larger than the SFR derived from the 
optical. Based on the sensitivity of our 24 /im data we 
infer an upper limit of 10 11 L Q for its total IR luminosity. 

Note that the possible association of the GRB 970508 
host with an ultra-luminous infrared galaxy proposed by 
Hanlon et al. (2000) based on ISOPHOT observations is 
clearly ruled out by the Spitzer data. 

4.3. The host galaxy of GRB 990705 

The GRB 990705 occured within a large and optically- 
luminous star-forming spiral galaxy observed face-on at 
z = 0.8424 (Le Floc'h ct al. 2002). The absolute mag- 
nitude of this GRB host corresponds to a 2L 4 galaxy 
at z~l. While GRBs preferentially occur within sub- 
luminous and young sources, the luminosity and the 
morphology of this object indicates therefore that some 
bursts can also take place within luminous and more 
evolved systems. 

The galaxy is clearly detected at 4.5 ^m with 
IRAC and 24 /xm with ' MIPS. Based on the MIPS 
detection we estimate a total infrared luminosity 
Lir = 1.8^5 xlO n L (see Table 1). Hence, this source 
belongs to the category of the Luminous Infrared Galax- 
ies. Assuming that the totality of the infrared emission 
is powered by star formation, its luminosity corresponds 
to a star formation rate SFR{ r ~ 32 M© yr -1 , which 
is somewhat larger than the SFR derived from the ob- 
served UV continuum (SFR UV ~ 5-8 M , Le Floc'h et al. 
2002). These characteristics and the observed morphol- 
ogy of the source are actually very similar to those of 
the IR-luminous spirals that were detected in the ISO 
and Spitzer deep surveys (e.g., Flores ct al. 1999; Zheng 
et al. 2004; Bell et al. 2005; Melbourne et al. 2005) and 
that dominate the star formation at z~l (e.g., Chary 
& Elbaz 2001; Le Floc'h et al. 2005). Contrary to many 
other GRB-selected galaxies, the host of GRB 990705 is 
therefore a very good representative of the sources re- 
sponsible for the bulk of the star-forming activity in the 
Universe at such rcdshifts. 

As shown in Figure 7, the lack of data from the /, J, 
H or K bands results in a rather poor sampling of the 
SED in the rest-frame optical and near-IR. This prevents 
us from infering the exact nature of the emission that 
dominates the luminosity of this galaxy at these wave- 
lengths. However the non-detection at 8/zm provides 
an additional strong constraint that firmly excludes the 
presence of hot dust dominating at short mid-IR wave- 
lengths. As illustrated by the comparison with the lo- 
cal templates, the IRAC detection at 4.5 /xm thus argues 
for an evolved and massive underlying stellar population 
dominating the near-IR emission. 



4.4. The host galaxy of GRB 970828 

The host of GRB 970828 appears as an early-stage in- 
teracting system at z = 0.96, with 3 components (ref- 
ered to as galaxies "A" , "B" and "C" ) located over a 
— 30kpc region (Djorgovski et al. 2001; Bloom et al. 
2001). The galaxy where the burst is believed to have 
occured (component "galaxy B" ) is faint at optical wave- 
lengths (R = 25.1) but it is one of the reddest GRB hosts 
(R — K = 3.6, Djorgovski et al. 2001) in the near-infrared 
sample studied by Le Floc'h et al. (2003). 

There is a clear detection of this merging system at 
4.5 /xm and 24itm. In the IRAC image it appears as a 
point source rough/ centered between the two galaxies 
"A" and "B" . These two components are only separated 
by -1.9" on the sky (Bloom et al. 2001).' They can 
not be resolved with Spitzer and it is not possible to 
unambiguously determine whether the infrared emission 
originates from only one or both of these objects. 

Based on the MIPS detection we estimate a total in- 
frared luminosity Lir = 1.4^QgXl0 11 LQ. This merging 
system is therefore a Luminous Infrared Galaxy like the 
host of GRB 990705, characterized by a dust-obscured 
star formation rate SFR ir ~ 24±g Mgyr" 1 . This SFR 
estimate is much larger than the one derived from the 
UV continuum or from the flux of the [Oil] emission line 
detected in the Keck spectra for the two galaxies "A" 
and "B" (SFR UV ~SFR [ou] -0.5-1 M yr-\ Djorgov- 
ski ct al. 2001). As shown in Fig. 7 the spectral energy 
distribution is steeply rising from the optical to the mid- 
infrared, which suggests the presence of a typical dust- 
enshrouded young starburst with no underlying evolved 
stellar population. 

Interestingly enough the GRB 970828 is generally ref- 
ered to as a typical "dark burst". The optical after- 
glow of this GRB was not detected despite a deep and 
prompt search down to — 24.5 mag and despite an ac- 
curate localization of the event thanks to the detection 
of its radio transient counterpart. Given the moderately- 
high redshift of its host galaxy, Djorgovski et al. (2001) 
thus suggested that its optical emission was likely sup- 
pressed by an intervening cloud of material within the 
host, invoking the dust extinction as one possible expla- 
nation for the origin of at least a fraction of these "dark 
GRBs" . Our infrared detection with Spitzer provides ob- 
viously a strong support for this hypothesis. In luminous- 
infrared galaxies, most of the UV radiation emitted by 
young and massive stars are indeed absorbed by dust and 
re-radiated at longer wavelengths. 

4.5. The host galaxy of GRB 980703 

This object is one of few GRB hosts character- 
ized by an unambiguous detection at radio wavelengths 
(Bcrgcr et al. 2001). Observations performed at the 
VLA after the fading of the GRB radio counterpart re- 
vealed at the location of the burst a persistent emis- 
sion of 68.0±6.6/xJy at 1.43 GHz. Assuming the far- 
infrared/radio correlation observed in local starburst 
galaxies (Condon 1992), Berger et al. (2001) have thus 
argued that the host of GRB 980703 is an ultra-luminous 
infrared galaxy (i.e., Lir — 10 12 L Q ) characterized by an 
SFR of several hundreds of M Q yr _1 . 

At the redshift of this object (z — 0.97, Djorgovski et al. 
1998), ULIRGs are very easily detected at the mid-IR 
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Spitzer wavelengths (i.e., IRAC + MIPS 24 /mi). Strik- 
ingly though, the GRB 980703 host galaxy is not detected 
in our data, except in the IRAC-4.5/mi channel. As 
shown in Fig. 7 its spectral energy distribution from the 
mid-IR to the radio is therefore not consistent with typi- 
cal starburst-dominated infrared-luminous galaxies, and 
our non-detection at 24 /im argues for a very modest star- 
formation rate (SFR- lr < 24 M Q yr _1 ) when compared to 
the SFR claimed by Berger et al. (2001). 

The excess of radio emission relative to the mid-IR 
flux can be quantified by the (724 parameter defined as 
(724 = logw{S2i^ml 'Syocm), where S\ is the observed 
monochromatic flux density at the considered wavelength 
(Appleton et al. 2004). Taking into account the sensi- 
tivity of our 24 /im data, we derive an upper limit of 
<724 = 0.09 for the host of GRB 980703, which is signifi- 
cantly lower than the average value (724 ~ 1 observed for 
starburst galaxies (Appleton et al. 2004). 

Low <724 measurements (i.e., (724 S^O) are usually inter- 
preted as the signature of radio-loud objects character- 
ized by an AGN-dominated SED (Higdon et al. 2005). 
Would the presence of an active galactic nucleus explain 
this radio detection in the host of GRB 980703 ? There 
is no evidence favoring this hypothesis in the optical 
spectrum of the host (Djorgovski et al. 1998), and no 
temporal variability in the radio emission has been re- 
ported by Berger et al. (2001). However we note that 
the slope of the radio continuum (spectral index /3~0.3, 
where F v oc v~P) is significantly flatter than the spec- 
trum characterizing the majority of the starburst " /iJy 
galaxies" (i.e., S^Ocm 100 /iJy ) selected at radio wave- 
lengths (/3~0.8, Richards 2000). As a result, the radio 
spectrum in itself would be more easily explained by the 
presence of an AGN than supernova remnants in star- 
bursting regions. In fact, the existing data do not al- 
low us to unambiguously disentangle between these two 
potential contributions, and this GRB host could still 
be characterized by a rare type of IR-luminous SED yet 
leading to a faint emission at 24 /im. Deep X-ray observa- 
tions and/or far-IR MIPS imaging at 70 /im and 160 /im 
should provide better constraints on the nature of this 
object and its spectral energy distribution. 

4.6. The host galaxy of GRB 980613 

The host of GRB 980613 is another merging system 
characterized by a very complex environment with up to 
9 galaxy fragments interacting with each other at z = 1.10 
(Chary et al. 2002; Hjorth et al. 2002; Djorgovski et al. 
2003). These components show a moderate activity of 
star formation in the optical (i.e., SFR opt . < 5 M© yr _1 ) 
but some of them display very red R — K colors. 

The system is clearly apparent in our Spitzer data. In 
addition to the detections at 4.5 /im and 24 /im it is actu- 
ally the only host galaxy of our sample also detected at 
8 /im. The flux estimated at 24 /im corresponds to a total 
infrared luminosity Lir = 5^g x 10 11 L , leading to an IR- 

cquivalent star formation rate SFR[ r ~ 87t-],2 6 M© y r_1 - 
However it should be noted that the Spitzer detec- 
tion does not coincide with the component of the inter- 
action where the GRB was observed (component "H", 
see Hjorth et al. 2002). It rather corresponds to the 
very red fragments denoted "C" and "D" by Chary et al. 
(2002) and Djorgovski et al. (2003). These two compo- 



nents are not spatially resolved by Spitzer but they are 
located more than 2.5" away from component "H" and 
they can be distinguished from the latter. We infer that 
the GRB 980613 did not occur in the region harboring 
the most intense star-forming activity of the system, as 
was already noted by Hjorth et al. (2002) based on deep 
optical HST data. 

The spectral energy distribution of the component de- 
tected with Spitzer in this interaction presents a striking 
contrast with the SEDs of the other GRB hosts observed 
in our data (see Figure 7). The fluxes measured at 4.5 /im 
and 8 /mi are particularly bright given the redshift of the 
host galaxy, and the 8 /mi detection reveals a clear in- 
flexion of the SED in the rest-frame near-infrared. This 
suggests not only a significant hot dust emission domi- 
nating the SED redward ~ 2 /im but also the contribution 
of an evolved and massive underlying stellar population 
likely dominating the optical wavelengths. 

4.7. The host galaxy of GRB 990506 

The GRB 990506 is another typical example of a "dark 
burst" undetected in the optical despite deep and prompt 
imaging after the gamma-ray explosion. Its accurate lo- 
calization was obtained based on the detection of its ra- 
dio afterglow (Taylor et al. 2000) and subsequent follow- 
ups revealed a host galaxy located at z = 1.31 (Bloom 
et al. 2003), characterized by a very compact morphol- 
ogy (Holland et al. 2000c) and a red R — K color 
(R - K~A, Le Floc'h et al. 2003). Its [Oil] emis- 
sion line indicates a dust-uncorrected star formation rate 
SFR [on] = 13 M yr" 1 (Bloom et al. 2003), which is sub- 
stantially higher than the median SFR characterizing 
the global population of the GRB host galaxies at optical 
wavelenghts. 

It is marginally detected at the 2a level at 4.5 /im, 
but it is not seen in the 8/iin and 24 /mi images. As 
already suggested by Barnard et al. (2003) who did 
not detect this object with SCUBA, this case reveals 
that "dark GRBs" are not systematically associated with 
dust-enshrouded massive star-forming activity. 

4.8. The host galaxy of GRB 010222 

The GRB 010222 has been widely refered as the pro- 
totypical burst associated with violent starburst activity 
at high redshift, bringing also further support for the 
death of young and massive stars as the origin of long 
GRBs. A persistent source with an average flux den- 
sity of 3.74 ± 0.53 mJy at 850 /mi and 1.05 ± 0.22 mJy at 
1.2 mm was indeed observed with SCUBA and MAMBO 
at the location of its afterglow (Frail et al. 2002). At 
the redshift of the burst (2; = 1.48, Jha et al. 2001), 
the reported flux densities correspond to an infrared lu- 
minosity Lir ~ 4x 10 12 Lq assuming typical far-infrared 
galaxy SEDs. Therefore Frail et al. (2002) argued that 
the GRB 010222 host is a dusty ultra-luminous infrared 
galaxy experiencing a very intense episode of star forma- 
tion (SFR~ 600 Moyr- 1 '). 

However, this object is not detected in our data (see 
Figure 8), which obviously raises some doubt regarding 
its association with the SCUBA/MAMBO source. The 
typical SCUBA and MAMBO galaxies have easily been 
detected with Spitzer (Charmandaris et al. 2004; Egami 
et al. 2004; Frayer et al. 2004; Ivison et al. 2004) while 
they are located on average at higher redshifts than the 
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host of GRB 010222 (e.g., Chapman et al. 2003). Fur- 
thermore, this GRB host galaxy is a very faint blue ob- 
ject that contrasts with the typical properties charac- 
terizing the optical counterparts of the SCUBA sources 
(Smail et al. 2004). Considering that it is not de- 
tected with more than 3er in any of the Spitzcr mid- 
IR bands nor at radio wavelenfgths (Berger et al. 2003), 
the claim for an association between the GRB host and 
the SCUBA/MAMBO source thus appears questionable. 
In fact Frail et al. (2002) mention the presence of 4 
other redder galaxies detected in the if-band and located 
within the 15"-diametcr beam of SCUBA centered at the 
position of the GRB host. Three of these galaxies are de- 
tected with IRAC and with MIPS, which could suggest 
that the SCUBA/MAMBO detection is likely associated 
with one or several of these other sources instead. 

Strictly speaking though, the non-detection of the host 
at 24 /im is not sufficient to completely rule out the pos- 
sible association between the GRB 010222 and a ULIRG. 
The SED of ultra-luminous infrared galaxies can be char- 
acterized by a strong silicate absorption at 9.7 /zm rest- 
frame (Spoon et al. 2004; Armus et al. 2004). At the 
distance of the host (z = 1.48), this feature would be red- 
shifted in the 24 /im band, which could explain the non- 
detection by MIPS at this wavelength (see Kasliwal et al. 
2005 for detailed discussion on this issue). In this case 
however, the non-detection at 4.5 /im would still remain 
very puzzling. 

4.9. The host galaxy of GRB 990123 

The high spatial resolution image of this galaxy ob- 
tained with the HST revealed a strongly-interacting sys- 
tem with a complex morphology, and the optical tran- 
sient of GRB 990123 was actually observed very near a 
star-forming region associated with one of its merging 
components (Bloom et al. 1999; Fruchter et al. 1999; 
Holland & Hjorth 1999). Located at 2 = 1.6 (Kulkarni 
et al. 1999), the GRB host has a luminosity of ~0.5L* 
in the optical and its observed UV continuum argues 
for a small amount of star formation (SFR ~ 4 M© yr -1 , 
Bloom et al. 1999). 

Our IRAC and MIPS observations did not lead to any 
detection in the infrared, which indicates that this sub- 
luminous galaxy is a low-mass quiescent starburst with 
no significant contribution from dust-enshrouded star- 
forming activity. 

4.10. The host galaxy of GRB 990510 

This host galaxy is a very faint object (U~28mag, 
Fruchter et al. 2000a) at z = 1.62 (Vreeswijk et al. 2001) 
with an absolute B-band magnitude Mb ~ -17. 20 mag 
(Le Floc'h et al. 2003). We did not detect it with IRAC 
nor with MIPS. This suggests that the faintness of this 
source observed in the optical does not originate from a 
high amount of dust extinction but rather points to an 
intrinsically low-mass and young object with negligible 
amount of star formation. 

4.11. The host galaxy of GRB 971214 

The host of GRB 971214 is one of the most distant 
objects that have been spectroscopically identified based 
on the optical transient of a GRB. Located at a redshift 
z = 3.418, it is an ~ L± galaxy with a somewhat irregu- 



lar morphology and a surface brightness probably dom- 
inated by an exponential profile (Kulkarni et al. 1998; 
Odewahn et al. 1998). Neglecting a possible extinction 
by dust, its UV continuum and Lya emission line both 
argue for a rather small star formation rate (SFR ~ 1 
SMqyt- 1 , Kulkarni et al. 1998). 

We did not detect this galaxy with Spitzer. At the 
redshift of the GRB, our 24 /zm data are only sensitive 
to Hyper Luminous Infrared Galaxies (Lir > 10 13 L Q ) at 
the flux limit of our observations. Therefore the non- 
detection by MIPS does not bring a strong constraint 
on the infrared properties of the host. At similar red- 
shifts however, typical Lyman-Break Galaxies (LBGs) 
are easily detected at the short IRAC wavelengths and 
some of them are also detected at 24 /tm (Barmby et al. 
2004; Huang et al. 2005). Hence, this lack of detection 
in our data indicates that the host of GRB 971214 is 
likely not as massive and evolved as the LBGs coexisting 
at the same epoch. It is not one of the most actively- 
starbursting LBGs neither. 

4.12. The host galaxies of GRB 980326, GRB 980329, 
GRB 980519 and GRB 990308 

The redshifts of these four GRB host galaxies have 
not been determined spectroscopically. Deep observa- 
tions performed with the HST have shown that they 
are very faint in the optical (i?>26mag, Jaunsen et al. 
2003; Fruchter et al. 2001), and none of them is actu- 
ally detected with Spitzer. This suggests that their faint 
optical emission is probably not due to a large extinc- 
tion by dust, but rather indicates galaxies with very low 
bolometric luminosities. In favor of this interpretation, 
we note that the very faint host galaxy of GRB 990510 
(V ^28 mag, z = 1.62) was neither detected in our sam- 
ple in spite of its redshift easily accessible for Spitzer 
(see Sect. 4.10). The possibility that these objects are at 
very high redshift can not be excluded though (see e.g., 
Fruchter 1999), in which case the non-detections with 
IRAC and MIPS would not be surprising. 

4.13. The host galaxy of GRB 981226 

As in the case of the four sources previously discussed, 
the spectroscopic redshift of the GRB 981226 host has 
not been established. However it is brighter in the op- 
tical (i?~24.5, Frail et al. 1999; Holland et al. 2000b; 
Saracco et al. 2001). It has been detected in the K s band 
(K s = 21.1 ±0.2, Lc Floc'h et al. 2003), and it is also 
clearly apparent in the IRAC 4.5 /xm image. Its R — K 
color makes it one of the reddest GRB host galaxies that 
have been studied so far, and its SED also looks rather 
steep between the K s band and the IRAC 4.5 /im chan- 
nel. 

Our current data set does not allow us to distinguish 
whether the IRAC detection and the red colors of this 
host reveals an old star population, a hot dust emis- 
sion, or a combination of both. Invoking the presence 
of dust may be attractive since the GRB 981226 is also 
one of the few "dark bursts" that did not exhibit any de- 
tectable optical counterpart. However we consider this 
explanation unlikely as this object is not detected at 
8.0 /im or 24 /im. Note that it was neither detected in 
the submillimctcr (i^so^m = 2.79 ± 1.17 mjy, Barnard 
et al. 2003) and the flux reported by Berger et al. (2003) 
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at radio wavelengths points to less than a 2a detection 
(-F8.46GHz = 21 ± 12 /xJy). The interpretation is further 
complicated by the fact that no spectroscopic redshift 
has ever been determined for this galaxy. 

5. DISCUSSION 

5.1. The origin of the infrared emission in the GRB 
host galaxies 

The IRAC and MIPS instruments on-board Spitzer 
have recently opened new exciting perspectives for tack- 
ling the evolution of galaxies in the early Universe. With 
unprecedented sensitivity at 3.6-8.0 ^m, the IRAC chan- 
nels can sample the rest-frame near-IR emission of very 
distant sources, while the MIPS imager can detect the 
hot dust emission of active starbursts and AGNs up to 
z~2-3. Both cameras provide therefore new insights to 
study the evolved stellar populations and the importance 
of dust-obscured star formation at high redshift. 

Interpreting the observed mid-IR emission of GRB- 
selected galaxies can be however more subtle than ana- 
lyzing the other mid-IR sources detected in the field. As 
mentioned in Sect. 3.2 Venemans & Blain (2001) have ar- 
gued that optical and UV afterglows can be significantly 
absorbed by dusty material surrounding GRBs, and that 
the reprocessed IR light can be observed several years 
after a burst because of the very slow time-scale varia- 
tions of the heated dust emission. The SED of this dust 
component would be characterized by a steep increase 
from the short wavelengths up to ~ 8-10 /zm rest-frame, 
followed by a gentle decline in the mid-IR and far-IR 
(see Fig. 5 of Venemans & Blain 2001). Up to redshift 
z~l-2, it could dominate the total luminosity of the 
host galaxy, then questioning the IRAC and MIPS flux 
measurements as mass and star formation rate indica- 
tors. 

However we did not find any obvious evidence for such 
"GRB-hcatcd" dust emission among the Spitzer-detected 
GRB hosts. In the resolved GRB 980425 host galaxy 
at z = 0.0085 we do not see any signature of this effect 
at the location of the hypernova SN 1998bw, while the 
infrared emission detected in the complex environment 
of the host of GRB 980613 (z = 1.10) docs not coincide 
with the galaxy fragment where the burst occured. At 
the redshift of the other sources (i.e., 0.84 < z <1.1), the 
IRAC-8/xm and MIPS-24^m bands constrain the rest- 
frame 4/zm/12^m flux ratios, which appear to be too 
red compared to the SED of the GRB dust component 
predicted by the models. 

The detectability of this GRB dust emission obvi- 
ously depends on the efficiency of the absorption of the 
UV/optical afterglow. As already noted by Venemans & 
Blain (2001), it should be preferentially observed toward 
the host galaxies of those "dark bursts" originating from 
dust-enshrouded star-forming regions 12 . Even though 
three of our galaxies fall indeed in this category (i.e., 
the hosts of GRB 970828, GRB 981226 and GRB 990506), 
our sample has been mostly built from GRBs pinpointed 
with optical afterglows, and our selection could also in- 
duce a bias against the detection of this burst-remnant 

12 The so-called "dark bursts" can also be GRBs with intrin- 
sically faint afterglows (e.g., Fox et al. 2003) or very high red- 
shift bursts with optical counterparts suppressed by hya absorp- 
tion (e.g., Lamb & Reichart 2000). 



dust emission. 

As a result, we conclude that the contribution of this 
component is negligible in our data. The IRAC and 
MIPS detections should truly reflect the properties of the 
stellar populations and the global star-forming activity 
within the observed host galaxies. 

5.2. A panchromatic view on the nature of the GRB 

hosts 

A few sources from our sample display clear sig- 
natures of evolved stellar populations and intense 
starbursting activity. An example is the host of 
GRB 990705 at z = 0.84, which harbors active star for- 
mation (SFRi T ~ 32 Mq yr _1 , see Sect. 4.3) and appears 
as a large and massive Sc spiral galaxy typical of the 
disk-dominated systems at these redshifts. Other cases 
of active starbursts, like the host of GRB 000418, have 
previously been found at radio and submillimeter wave- 
lengths (Berger et al. 2003). 

On average though, our MIPS observations combined 
with the submillimeter and radio published photometry 
do not really favor a population of GRB host galax- 
ies dominated by very active and luminous dusty star- 
bursts. With the exception of the host of GRB 980703 
that we discussed in Sect. 4.5, our 24 /xm non-detections 
are consistent with the existing SCUBA and VLA data 
and they are even more constraining if one assumes typ- 
ical IR-luminous and starburst-dominated SEDs for the 
host of GRBs (see Figure 7). Our measurements argue 
against the presence of numerous LIRGs and ULIRGs in 
the GRB host population and they also point to lower 
star formation rates (see Sect. 3.2 and Table 1) which 
are not consistent with some conclusions obtained by 
other groups. For instance Berger et al. (2003) have 
recently claimed that 20% of GRB host galaxies have 
SFR ~ 500 M Q yr _1 , which is obviously in disagreement 
with our Spitzer results. 

Similarly our IRAC data argue very clearly against a 
population of sources with massive and evolved stellar 
populations dominating their near-IR rest-frame emis- 
sion. Assuming a standard conversion between the mass 
and the rest-frame near-IR absolute luminosity (i.e., - 
0.25 < logm (M/L K ) < +0.15, Bell et al. 2003), the 
fluxes or upper limits that we measure at 4.5 /xm translate 
for most cases into masses M. < 5xl0 9 M Q . The IRAC 
data are also consistent with the constraints derived from 
the optical and the near-IR published photometry as- 
suming typical SEDs of blue star-forming galaxies (see 
Figure 7), and our results globally agree with the rather 
small masses determined by Chary et al. (2002) using 
deep if-band observations of GRB hosts at Keck. Our 
interpretation, on the other hand, strikingly contrasts 
with the conclusions recently obtained by Conselice et al. 
(2005) who argue for a trend toward massive sources at 
z > 1 based on an analysis of the GRB host morphologies. 

Hence the Spitzer view on the GRB-selccted galaxies 
rather suggests low mass sources characterized by a rel- 
atively modest amount of dust obscuration. Similar con- 
clusions have already been derived based on the optical 
and near-IR properties of these objects (e.g., Le Floc'h 
et al. 2003; Courty et al. 2004; Christensen et al. 2004). 

5.3. Implications to GRBs as SFR probes 
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Because of their dust-penetrating power and their re- 
lation with the death of young and massive stars it has 
been argued that the long GRBs could be used as an 
unbiased probe of the star formation history of the Uni- 
verse. On a pure statistical point of view, their host 
galaxies should be therefore representative of the sources 
producing the bulk of the stellar mass at high redshift. 
Arc our results consistent with this picture ? 

The detection of the cosmic infrared back- 
ground and the recent surveys performed at in- 
frared/submillimeter/radio wavelengths have revealed 
that a significant fraction of the present-day stellar 
mass was formed in the past during short-live and 
dust-obscured episodes of intense star formation within 
infrared-luminous galaxies (e.g., Puget et al. 1996; 
Blain et al. 1999; Chary & Elbaz 2001; Cowie et al. 

2004) . At z~l, these infrared- luminous starbursts 
(i.e., Lir>10 11 L q ) are responsible for ~70% of the 
star-forming activity in the Universe (Le Floc'h et al. 

2005) and their contribution is believed to be even more 
important at higher rcdshifts (e.g., Blain et al. 2002; 
Lagache et al. 2003). Using Spitzer-upd&ted IR galaxy 
evolution models (e.g., Lagache et al. 2004; Chary 
et al. 2004) we estimate that more than half of the 
global star-forming activity throughout the lifetime of 
the Universe has occured within galaxies characterized 
by /24um 100 /iJy In this context, the significant 
fraction of non-detections in our MIPS data is therefore 
surprising. 

Similarly, the downsizing evolution of the cosmic star 
formation history reveals that the bulk of the star- 
forming activity has been moving from massive galaxies 
at high redhifts to low-mass objects in the present-day 
Universe (e.g., Cowie et al. 1996; Juneau et al. 2005). 
Should the long GRBs trace the whole population of 
distant starbursting sources, we should therefore detect 
these bursts preferentially toward massive systems easily 
accessible to IRAC. 

This apparent discrepancy between the nature of the 
GRB host galaxies as a whole and the sources dom- 
inating the high-redshift star-forming activity has al- 
ready been noted from a comparison between their col- 
ors and luminosities at optical and near-IR wavelengths 
(e.g., Le Floc'h et al. 2003) or their star formation rate 
recovered from cosmological simulations (Courty et al. 
2004). For instance GRB hosts appear blue and optically 
sub-luminous compared to the massive dusty starbursts. 
Very often they also display US in emission (Djorgov- 
ski et al. 1998; Soderberg et al. 2004; Prochaska et al. 
2004; Gorosabel et al. 2005), which reveals star forma- 
tion episodes characterized by younger populations than 
those typically observed in the distant luminous-infrared 
galaxies (Flores et al. 1999; Hammer et al. 2001, 2005). 
As a result we conclude that the GRB afterglows, as they 
are currently selected, can not be considered as unbiased 
probes of the integrated activity of star formation in the 
high redshift Universe. 

5.4. A bias in the Spitzer sample ? 

As we have described in more detail earlier, our sam- 
ple consists mostly of galaxies pinpointed using optical 
afterglows identified prior to 2001, when GRB follow-ups 
were not as prompt and efficient as they are currently 
This could bias the selection toward the most luminous 



transients, and therefore against dust-enshrouded GRBs 
occuring within luminous-infrared galaxies. A relevant 
case illustrating this hypothesis is given by the dark burst 
GRB 970828 and its host galaxy (see Sect. 4.4). The ac- 
curate coordinates of this GRB were determined thanks 
to the detection of its radio afterglow, and the absence 
of optical transient emission was interpreted as an evi- 
dence for dust extinction in the host galaxy (Djorgov- 
ski et al. 2001). The detection of this object at 24 /im 
clearly supports this interpretation, and it suggests that 
a sample of GRB hosts purely selected with optical af- 
terglows could be biased toward dust-free sources. In 
this context, the XRT instrument on-board the recently- 
launched Swift satellite is now routinely localizing X-ray 
GRB afterglows with an accuracy of ~ 5" on the sky. In 
the near-future, the comparison between the properties 
of GRB host galaxies selected from optical and X-ray 
afterglows will provide better insights into this possible 
bias affecting our current sample. 

On the other hand, we have not detected the two 
other host galaxies selected from radio afterglows with 
no optical transient (i.e., the hosts of GRB 981226 and 
GRB 990506), which shows that not all dark GRBs orig- 
inate from dusty galaxies. A similar conclusion was 
reached by Barnard et al. (2003) based on the non- 
detection of four dark GRB hosts with the SCUBA cam- 
era at 850 /im. In fact, many of these dark bursts could 
simply be associated with intrinsically very faint and/or 
fast-decaying optical afterglows (e.g., Fynbo et al. 2001; 
Fox et al. 2003). 

In addition to this observational bias that could af- 
fect our selection, another explanation might be directly 
related to the physical properties characterizing the lo- 
cal environments where long GRBs take place. Based 
on theoretical simulations and the connection between 
GRBs, hypernovae and Type Ic supernovae, it has been 
argued that long GRBs should more likely occur within 
binary systems (Izzard et al. 2004; Podsiadlowski et al. 
2004; Mirabel 2004a,b), whose frequency in star- forming 
regions may vary with redshift. Furthermore it has been 
suggested that long gamma-ray bursts would be more ef- 
ficiently produced - and they would also appear more 
luminous - if they originate from low metallicity progen- 
itors (e.g., MacFadyen & Woosley 1999; Ramirez- Ruiz 
et al. 2002; Meynet & Maeder 2005; Hirschi et al. 2005). 
Massive stars with metal-poor envelopes keep a high an- 
gular momentum in the latest stage of their evolution, 
and they are also less subject to mass loss. After the 
final collapse, this would favor the formation of a fast- 
rotating black hole with accretion of material, that could 
more easily lead to a bright GRB event. In this case 
GRBs would be preferentially observed in young and 
chemically-unevolved galaxies, which would explain this 
lack of massive and dusty starbursts in our Spitzer data. 

In fact, direct evidence for low metallicity has already 
been observed in several GRB host galaxies (Soderberg 
et al. 2004; Prochaska et al. 2004; Gorosabel et al. 2005). 
Furthermore, the global properties of these objects such 
as their blue colors, their relatively low luminosities, their 
Lya emission and their high specific star-formation rate 
(Le Floc'h et al. 2003; Fynbo et al. 2003; Christensen 
et al. 2004) do support a picture where GRBs occur 
in low-mass, young and hence metal-poor starbursts. 
This is also corroborated by the cosmological Simula- 
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tions obtained by Courty et al. (2004) who identify the 
GRB hosts as the most efficient star- forming objects (i.e., 
sources with the highest specific star-formation rate) and 
not as galaxies with obvious high SFR. In fact a good 
illustration of this global property can be given by the 
complex environment of the host of GRB 980613. This 
burst did not occur in the region detected by MIPS and 
harboring therefore the most intense star-forming activ- 
ity within the galaxy (see also Hjorth et al. 2002 for a 
similar conclusion based on optical HST data). Note that 
the same interpretation can be derived from the charac- 
teristics of the GRB 980425 host galaxy, since the hyper- 
nova SN 1998bw did not occur in the most active region 
as revealed by the very luminous mid-IR point source 
detected with Spitzer. 

6. SUMMARY 

We have presented 4.5^m/8.0^m-IRAC and 24 ^m- 
MIPS observations of 16 GRB host galaxies, and our 
results can be summarized as follows: 

• It is now well-established that long GRBs arc mark- 
ers of recent bursts of star formation in galaxies. 
However most of the GRB hosts in our sample were 
not detected in the rest-frame near-IR and mid- 
IR with Spitzer, which argues against a popula- 
tion of sources globally dominated by massive and 
IR- luminous starbursts (Lir > 5xlO n L Q ). Cur- 
rent IR galaxy evolution models indicate that more 
than half of the integrated star-forming activity 
throughout the lifetime of the Universe occured 
within /24/jm <; 100 /iJy galaxies easily detectable 
by MIPS-24/im. In this context our results imply 
that GRBs identified and localized with the current 
techniques can not be used as unbiased probes of 
the global star formation in the early Universe. 

• The detection of the GRB 970828 host by MIPS 
at 24 fxm strongly supports the idea that some of 
the so-called "dark bursts" can be explained by the 
effect of dust extinction within their host galaxy. 
Even though the hosts of the two other dark GRBs 
in our sample (GRB 981226 and GRB 990506) were 
not detected, this could indicate that the currently 
GRB-selected sources are biased against dusty star- 
bursts. In the near future the localization of GRB 
hosts using X-ray afterglows detected with the 
XRT instrument on-board Swift will provide new 
insights on this potential bias. 
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• The host of GRB 010222 that has been claimed 
to be associated with a SCUBA/MAMBO galaxy 
at z = 1.48 is not detected with IRAC nor with 
MIPS, thus bringing strong doubt on the identi- 
fication of this host at long wavelengths. Similarly 
our non-detection of the host of GRB 970508 at 
z = 0.84 rules out the proposed association of this 
object with an ultra-luminous infrared galaxy seen 
at 90 /xm with ISOPHOT, and the non-detection 
of the host of GRB 980703 at 24 ^m suggests that 
the radio emission previously-detected in this ob- 
ject does not originate from massive star formation. 

• The observations of the host galaxies of 
GRB 980425 and GRB 980613 reveal that these 
bursts did not occur in the regions harboring 
the most active star-forming activity within 
their hosts. This favors a picture in which the 
production of GRBs does not exactly scale with 
star formation but also depends on other pa- 
rameters that remain to be explored (e.g., age 
and chemical enrichment of the parent progenitor 
populations, Initial Mass Function, fraction of 
binary systems,...). 
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TABLE 1 

MlD-IR PHOTOMETRY OF GRB HOST GALAXIES 



Name 


RA(J2000) 


Dec(J2000) 


(**Jy) 


f 8 „™ (a) 

(/*Jy) 


(/*Jy) 


Redshift 
z Reference 


Lir CO 
(L©) 


(M@ yr- 1 ) 




06: 


53 


:49.45 


4-7Q-1 Q ^ 


< . o.u 


±o.o 


' OO 


0.84 


1 


<r i y i nil 


-\ 7 
v. ± ( 




18: 


:08 


31.60 


+oy:18:ol.o 


3.7+0.2 


< 17.8 


85+15 


0.96 


2 


1.4_ g x 10 


Ozl+43 
24_ 14 


GRB971214 


11: 


56 


:26.40 


+65:12:00.5 


< 3.6 


< 17.1 


< 92 


3.42 


3 


< 4.7 X 10 


< 8100 


vjitr>9oUozo 


08: 


:36 


:34.28 


-lo.oi.zo.y 


^ Q P. 


/ it; 1 
< 10.1 


< OD 










GRB980329 


07: 


:02 


38.02 


+38:50:44.0 


< 4.8 


< 24.4 


< 97 










GRB980425 


19: 


35 


:03.02 


-52:50:44.8 


2 950+100 


11900+300 


27300+200 


0.0085 


4 


2 x 10 9 ( d > 


0.4 


GRB980519 


23: 


:22: 


:21.50 


+77:15:43.3 


< 4.2 


< 24.4 


< 93 










GRB980613 


10: 


:17 


:57.82 


+71:27:25.5 


36.1+1.6 


34.9+1.7 


170+30 


1.10 


5 


5±f x 10 11 


87±^ 6 


GRB980703 


23: 


:59: 


:06.67 


+08:35:07.1 


11.2+0.6 


< 23.7 


< 83 


0.97 


6 


< 1.4 x 10 11 


< 24 


GRB981226 


23 


29 


:37.21 


-23:55:53.8 


4.1+0.2 


< 29.6 


< 90 










GRB990123 


15: 


:25 


: 30.31 


+44:45:59.2 


< 3.6 


< 17.1 


< 82 


1.60 


7 


< 8 x 10 11 


< 140 


GRB990308 


12: 


:23: 


:11.44 


+06:44:05.1 


< 4.8 


< 23.7 


< 87 










GRB990506 


11: 


:54: 


:50,14 


-26:40:35.0 


1.6 (e) 


< 23.7 


< 80 


1.31 


8 


< 4.8 x 10 11 


< 83 


GRB990510 


13: 


:38: 


:07.11 


-80:29:48.2 


< 4.2 


< 18.4 


< 98 


1.62 


9 


< 9.8 x 10 11 


< 170 


GRB990705 


05: 


:09: 


:54.50 


-72:07:53.0 


17.2+0.8 


< 17.8 


150+20 


0.84 


10 


1.81q6 x 1Ql1 


32t^ 


GRB010222 


14: 


:52 


: 12.55 


+43:01:06.3 


< 3.0 


< 21.1 


< 81 


1.48 


11 


< 7.4 x 10 11 


< 130 



All upper limits indicated arc 3c. 
(k) defined as the energy density integrated between 8 fxra and 1000 fim, and computed assuming a ACDM cosmology with 

H =70 km s" 1 Mpc" 1 , Q m =0.3 and Q x = 0.7. 
( c ) assuming SFR JR (M Q yr" 1 ) = 1.72 x 10~ 10 L IR (L ). 

More than 75% of the 24 /xm monochromatic luminosity of this galaxy originates from a single HII region (see Sect. 4.1). This unusual 
property makes the extrapolation to the total infrared highly uncertain. 

( C J The source is only detected at a ~ 2c level. 

References: (1) Bloom ct al. (1998) ; (2) Djorgovski ct al. (2001) ; (3) Kulkarni et al. (1998) ; (4) Tinney et al. (1998) ; (5) Djorgovski et al. (2003) ; (6) 
Djorgovski et al. (1998) ; (7) Kulkarni ct al. (1999) ; (8) Bloom ct al. (2003) ; (9) Vrccswijk ct al. (2001) ; (10) Lc Floc'h ct al. (2002) ; (11) Jha ct al. (2001). 
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Fig. 1 — 
each image, 
to the left. 



Postage stamp mosaic of the area surrounding the GRB 970828 host (z = 0.96). The location of the GRB is at the center of 
a) Image at 4.5/j.m, b) Image at 8.0/im, c) Optical HST image (see text), d) Image at 24fim. North is to the top and East is 
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Fig. 2. — Postage stamp mosaic of the area surrounding the GRB 980425 host at 2 = 0.0085 (same legend as in Fig. 1). Note that the 
bright point source detected at the Spitzer wavelengths is not exactly located at the position of the GRB but rather coincides with the 
HII region detected 6" away in the north-west direction. 
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Fig. 3. — Postage stamp mosaic of the area surrounding the GRB 980613 host at 2 = 1.10 (same legend as in Fig. 1). 
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Fig. 4. — Postage stamp mosaic of the area surrounding the GRB 980703 host at z = 0.97 (same legend as in Fig. 1). 
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Fig. 6. — Postage stamp mosaic of the area surrounding the GRB 990705 host at z = 0.84 (same legend as in Fig. 1). 
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Fig. 7b. — Monochromatic luminosities of GRB host galaxies computed from the optical to the radio (red filled diamonds, upper limits 
appear as red downward arrows). Optical, near-IR, submillimeter and radio photometry measurements are taken from Bloom ct al. (1998), 
Kulkarni ct al. (1998), Bloom et al. (1999), Fruchter et al. (1999), Holland & Hjorth (1999), Vreeswijk et al. (1999), Fruchter et al. (2000b), 
Fynbo et al. (2000), Berger et al. (2001), Djorgovski et al. (2001), Holland et al. (2001), Sokolov et al. (2001), Chary et al. (2002), Frail 
et al. (2002), Le Floc'h et al. (2002), Barnard et al. (2003), Berger et al. (2003), Djorgovski et al. (2003), Le Floc'h et al. (2003) and 
Tanvir et al. (2004). For comparison we also show the SEDs of NGC 253 (orange), M 82 (blue) and Mrk 231 (black) as well as the SED of a 
starburst-dominated ULIRG with Ltr = 4 X 10 12 L Q (light brown, denoted "SF-ULIRG") and the SED of a cold LIRG with Li R = 10 11 L Q 
(purple). All these SEDs were corrected for the distance and they are displayed in units of luminosity. As indicated in each panel, a scaling 
factor was applied to some of them in order to match the observed GRB host photometry at certain wavelengths. The bands where this 
normalization was performed were chosen on a case-by-case basis to better highlight the differences in luminosity at the other wavelengths 
(see text for more detail). Note that the host of GRB 980613 is a merging system, and the photometry reported in the figure does not refer 
to the region where the GRB occured but only to the component detected by Spitzer (see Sect. 4.6). 
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Fig. 8. — Postage stamp mosaic of the area surrounding the GRB 010222 host galaxy at 2 = 1.48 (same legend as in Fig. 1). The 
15"-diameter beam of SCUBA centered on the host is displayed in the IRAC 4.5 /im and MIPS 24mic images. In panel a), we also indicate 
the three sources refered to as galaxies "A", "B" and "C" by Frail et al. (2002). 



